function xk=cyclic41(s_in,x_k1)

byte_sequence=qpsk2bin(s_in);
g=[exp(-pi*j/2) 0; 0 exp(-pi*(j/2))];

gk=g^(byte_sequence(2)*2 + byte_sequence(1));
xk=gk*x_k1;
